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ABSTRACT 

We present a comparative study of the thermal emission of the transiting exoplanets WASP- lb and WASP- 
2b using the Spitzer Space Telescope. The two planets have very similar masses but suffer different levels of 
irradiation and are predicted to fall either side of a sharp transition between planets with and without hot strato- 
spheres. WASP- lb is one of the most highly irradiated planets studied to date. We measure planet/star contrast 
ratios in all four of the IRAC bands for both planets (3.6-8.0 ^m), and our results indicate the presence of a 
strong temperature inversion in the atmosphere of WASP- lb, particularly apparent at 8 /im, and no inversion in 
WASP-2b. In both cases the measured eclipse depths favor models in which incident energy is not redistributed 
efficiently from the day side to the night side of the planet. We fit the Spitzer light curves simultaneously with 
the best available radial velocity curves and transit photometry in order to provide updated measurements of 
system parameters. We do not find significant eccentricity in the orbit of either planet, suggesting that the 
inflated radius of WASP- lb is unlikely to be the result of tidal heating. Finally, by plotting ratios of secondary 
eclipse depths at 8 /im and 4.5 /im against irradiation for all available planets, we find evidence for a sharp 
transition in the emission spectra of hot Jupiters at an uradiation level of 2 x lO^ergs"' cm"^. We suggest 
this transition may be due to the presence of TiO in the upper atmospheres of the most strongly irradiated hot 
Jupiters. 

Subject headings: 



1. INTRODUCTION 

The Spitzer Space Telescope dWerner et al.l 12004 *) has 
been used to carry out the first photometry and emis- 
sion spectro scopy of exoplane t s that orbit main se- 
quence stars (iDeming et al.l 120051: iCharbonneau et al.l l2005t 
lRichardsonetal.1 120071: iGriUmair et all l2007h . This was 
achieved by observing transiting planets at secondary eclipse 
(when the planet is eclipsed by the star) which allows the 
emission of the planet to be separated from that of the 
star Photometry and spectroscopy are the key measurements 
needed to determine the physical properties of any astronom- 
ical object, and secondary eclipse observations allow us to 
consider the temperature structure and chemical composi- 
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tion of exoplanet atmospheres and the redistribution of en- 
ergy from t he day si de to the night side of the planet (e.g . 
Burrows et a l. 2005; Fortnev et all 120051: ISeager et all 120051: 
Barman et al.ll2005i) . 

Secondary ecUpse detec tions have been made 
from the ground (e.g. Ide Mooij & SnellenI 120091: 
iSing & Lopez-MoralesI 120091) but since the signal is 
weak even in the best cases (of order 0.1 per cent in hot 
Jupiters) the bulk of measurements to date are from space. 
Togeth er with detections in the optical with CoRoT and 
Keple r dSnellen et al.ll2009l |Alonso et al."2009; Borucki et alj 
120091) and the near infrared with i/5r (Swain et al, 2009a,bi), 
Spitzer is providing an increasingly clear picture of the 
thermal emission of exoplanets. 

Spitzer detections of thermal emission have been reported 
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Target Prog. 



WASP-lb 30129 
282 

WASP-2b 30129 
282 



Date 



2007-09-08 

2006- 12-30 

2007- 07-01 
2006-11-28 



Table 1 

Log of Spitzer observations of WASP-lb and WASP-2b. 



Start time 
UTC 

14:02:56 
14:15:27 
13:44:18 
08:23:11 



Duration Bary. corr. 



27 586 
27 146 
12427 
12032 



-11.8 



No. of frames X effective exposure per frame 
3.6 ^m 4.5 /im 5.8 /im 



1.0 iMm 



910 X 10.4 s 



Pipeline 
version 



+238.2 2072X 10.4s 2072x 10.4s S16.1.0 

+252.5 2055 X 10.4s 2055 x 10.4s S15.0.5 

+240.3 1818 X 1.2s 909 x 10.4s S16.1.0 



910 X 10.4s S15.0.5 



in various combinations of wavebands between 3.6 and 
24 /im. The bulk of the newly-discovered hot Jupiters are 
detectable only in the shorter wavelength S pitzer hands of the 
IRAC instrument (3.6, 4.5, 5.8 & 8.0 fim: iFazio et al.ll2004 
but fortunately this is a range in which strong molecular bands 
are expected, providing good constraints on atmospheric con- 
ditions. The IRAC observing modes also allow more efficient 
observations of these fainter systems, to some extent compen- 
sating for their lower brightness. The rapidly growing number 
of moderately-bright transiting hot Jupiters therefore provide 
an excellent opportunity to improve our understanding of ex- 
oplanet atmospheres. 

The existing Spitzer observations show that hot Jupiters 
are strongly heated by their parent stars, with typical bright- 
ness temperatures in the range 1 000-2000 K, and that their 
spectra deviate strongly from black bodies. Preliminary the- 
oretical calculations predicted strong molecular absorption 
in the IRAC bands feurr ows et al. 2005; Fortney et al. 200l 
ISeager et al.ll2005l:lBarman et a l. 2005), for which ther e is ev- 
idenc e in some systems (e.g. HD 189733b; Charbonneau^LalJ 
120081) . however other systems have measured brightness tem- 
peratures in excess of expectations, indicating emission fea- 
tures in the IRAC bands (e.g. HD209458b; iKnutson et aP 
l2008h . In these cases it is thought that an opacity source 
high in the atmosphere results in a hot stratosphere and that 
the molecular ban ds are driven into emission by the tem- 
perature inversion (iHube nv e t all 120031; iFortnev et al.l l2006t 
Harrington et al. 2007; Burro ws et al.ll20 07; Sing et al. 2008). 
Fortney et al.. (i2008i) and iBurrows et alJ(l2008 ) suggest that ir- 
radiated exoplanets may fall into two distinct classes, those 
with and those without hot stratospheres, depending on 
the level of incident ste llar flux ( dubbed pM and pL class 
planets respectively by IFortnev et al. 2008). The bright- 
est and best studied systems, HD209458b and HD 189733b, 
fall eith er side of the pred icted transition between these 
classes (IFortnev et al.l l2008h . and their IRAC fluxes sup- 
port the prese nce o f a ho t stratosphere in HD209458b 
jBurrows et aljp 007; Knutson et al. 2008) and its absence in 
HD 189733b dCharbonneau et al. 2008) . Most of the other 
systems that have been observed in a ll four IRAC bands 
also seem to support this overall picture (Knutson et al."2009t 
Machalek et al. 2009; Todorov et al. 2010; O' Donovan et all 
20101; iMachalek et al.ll2010t |Campo et al. 2010), but XO-lb 
and TrES-3 do not, with XO-lb presenting evi dence for a 
tempe rature inversion despite low irradiation ( Ma chalek et al.l 
120081 12009) . and TrES-3 not exhibiting e vidence for a temper - 
ature inversion despite high irradiation (iFressin et al.ll2010l) . 
It may be that additional parameters dictate the presence of a 
temperature inversion, although to some extent the picture is 
confused by the different models an d criteria for i nversi on de- 
tection applied by different authors. iGillon et"an (120101) show 
that the planets studied to date cover a wide region in color- 
color space, and do not fall clearly into two groups. 

In this paper we present Spitzer IRAC secondary eclipse 



detections of the ti-ansiting planets WASP-lb and WASP-2b, 
which were disc overed by the Wide Ang le Search for Planets 
(WAS P) project (IColHer Cameron et al.l2007b; Pollacco et alj 
HOO^. These planets make an interesting pair for com- 
parative study since they have near identical mass (0.9Mj) 
and yet WASP-lb is highly irradiated and expected t o have 
a hot stratosphere, while WASP-2b is not ( Fortne y" et al.l 
120081) . Indeed, WASP-lb is one of the most highly irra- 
diated planets studied with Spitzer to date (incident flux of 

2.5 X 10^ erg s~' cm"^). It is also one of the group of oversized 
hot Jupiters that h ave radii larger than can b e explained with 
canonical models dCharbonneau et an2007l) . 

In addition to presenting and discussing Spitzer detections 
of WASP-lb and WASP-2b, we present revised parameters 
for both planets based on simultaneous fits of all available 
photometry and spectroscopy. 

2. OBSERVATIONS 

2.1. Spitzer secondary-eclipse observations 

The Spitzer IRAC instrument provides images of two adja- 
cent fields, each in two w avebands (one fiel d at 3.6 & 5.8 fim, 
the other at 4.5 & 8.0 ^m;l Fazio et alj|2004l) . We observed the 
expected times of two secondary eclipses for each of WASP- 
lb and WASP-2b, allowing us to cover all four wavebands 
without repointing the telescope during an observation. A log 
of our four observations is given in Table[T] 

The precise target positions were carefully chosen in order 
to avoid bad pixels, keep saturated stars off of the array, ex- 
clude bright stars from regions known to scatter light onto the 
IRAC detectors, and to place bright comparison stars on the 
detector (although these were not used in our final analysis). 
The pointing position was not dithered during the observa- 
tions. 

Our observations were made in full array mode with 
12 s frame times (10.4 s effective exposure), except for the 
3.6/5.8 fim observation of WASP-2b, which was made in stel- 
lar photometry mode with pairs of 2 s frames taken in the 

3.6 /im band for each 12 s frame in the 5.8 /im band. This 
mode was used in order to avoid saturating the target in 
the 3.6 /tm band. Observation durations are listed in Ta- 
ble[T] and were chosen in order to cover approximately twice 
the expected secondary eclipse duration (transit durations 
are 3.7h and 1.8h for \V ASP-lb and WASP-2b respectively; 
ICharbonneau et al.ll2007l) . 

2.2. Radial velocity and transit observations 

In addition to the Spitzer secondary eclipse observations 
we also re-analyzed radial-velocity measurements of WASP- 
1 and WASP-2 from Coflier Cameron et al. (2007b) a nd z- 
band p rimary transit observations from ICharbonnea u et alJ 
(120071) as well as the SuperWA SP discovery photometry from 
IColher Cameron et all (l2007bl) . 

For WASP-1 we also included the /-band transit observa- 
tions of .Shporer et al.. (.2007.) in our analysis. 
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Figure 1. Upper panels show the Spitzer IRAC light curves of the star WASP-1 during the expected times of secondary eclipse of the planet WASP-lb. Points 
show the measurements from individual images, crosses show the data binned into five hundred b ins p e r orb ital period. Solid lines show our best fits to the 
echpse fight curves, including finear decorrelation and other instrumental effects described in Sects. I4n & l4. 21 The lower panels show the radial velocity curve 
of WASP-1 (measured with SOPHIE at OHP) and the z-band optical figfit curve during primary transit (from Keplercam), eacfi plotted on different pfiase ranges. 
These data were fitted simultaneously with the secondary eclipse observations. 



For WASP-2 we included additional radial-velocity mea- 
surements made with the Swiss 1 .2-m Euler telescope at La 
Silla, Chile between 2008 October 14 and November 6, and 
six pairs of measurements made with HARPS on the ESO 3.6- 
m telescope at La Silla on the nights of 2008 October 14/15, 
15/16 and 16/17. 

All the radial-velocity measurements used in our analy- 
sis are plotted in the lower-left panels of Figs.[T]&|2l Pri- 
mary transit photometry is shown in the lower-right panels 
of these figur es, but we only plot the transit photometry from 
ICharbonneau et al.i (120071) for clarity. 

3. 5pteer DATA REDUCTION 

For our analysis we used the Basic Calibrated Data (BCD) 
files generated by the IRAC pipelineQ One image is produced 
for each exposure, with the processing including: dark current 

' Pipeline description document: |http://ssc. spitzer caltech.edu/irac/dh/PDD.pdf] 



and bias subtraction; muxbleed correction for the InSb arrays; 
non-linearity correction; scattered light subtraction; flat field- 
ing; and photometric calibration. The pipeline version used to 
process each observation is given in Table[T] 

We first renormalized each BCD image into units of elec- 
trons by multiplying by exposure time and detector gain and 
dividing by the flux conversion factor given in the BCD file 
headers. 

Light curves were extracted using simple aperture photom- 
etry with a source aperture radius of four pixels in all bands. 
The aperture position was centred on the target in each image 
by centroiding the x and y pixel positions within a 9 pix search 
box. The background contribution was estimated by calcu- 
lating the mode of pixel values in an annulus centred on the 
source with inner and outer radii of 8 and 20 pixels. Uncer- 
tainties were estimated using counting statistics normalised to 
the measured variance in the background in each image. 

In order to remove contamination of the light curves by 
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Orbital phase Orbital phase 

Figure 2. Upper panels show the Spitzer IRAC light curves of the star WASP-2 during the expected times of secondary eclipse of the planet WASP-2b. Points 
show the measurements from individual images, crosses show the data binned into five hundred b ins p e r orb ital period. Solid lines show our best fits to the 
eclipse light curves, including linear decorrelation and other instrumental effects described in Sects. I4n & l4. 21 The lower panels show the radial velocity curve 
of WASP-2 (measured with SOPHIE, Coralie and HARPS) and the z-band optical light curve during primary transit (from Keplercam), each plotted on different 
phase ranges. These data were fitted simultaneously with the secondary eclipse observations. 



radiation hits we adopted a two stage process in which we 
first rejected images for which photometric measurements 
were extreme outhers from the mean of the entire dataset 
(> lOcr). We then formed a running mean of the remaining 
measurements across one hundred images and rejected im- 
ages that were highly significant outliers from the running 
mean (> 5a). In total, 31 1 of 12,806 images were rejected (2 
per cent). The proportion of frames rejected from each light 
curve is given in Table|2l 

The UTC times from the image headers were corrected to 
the Solar System barycenter using our own code and the co- 
ordinates of the Spitzer spacecraft from the HORIZONS on- 
line ephemeris system of the Jet Propulsion Laboratory^ The 
corrections applied to each observation are listed in Table[T] 
The maximum change in correction during an observation 
was 2 s, and so it was sufficient only to apply a single aver- 

^ http://ssd.jpl.nasa.gov/horizons.cgi 



age offset to each observation. 

The resulting Spitzer light curves are plotted in Figs.[T]&|2l 
The light curves were fitted at full resolution (described in the 
following section), but they are plotted here also binned into 
five hundred phase bins for clarity. 

4. LIGHT CURVE DECORRELATION AND PARAMETER FITTING 

4.1. 3.6 and 4.5 /^m decorrelation 

All four IRAC detectors exhibit distinctive patterns of cor- 
related systematic error. Data from the 3.6 /xm and 4.5 
/im (InSb) detectors are most strongly affected by intra- 
pixel variations in quantum efficiency (Reach et al.' '2005|; 
Charbonneau et al. 2005; Morales-Calderon et al. 2006). 
While the impact of these sensitivity differences is mini- 
mized by not moving the pointing during an observation, the 
spacecraft pointing tends to oscillate around the nominal posi- 
tion, with an amplitude of around 0. 1 pix, leading to position- 
dependent variations in the measured stellar flux. The effect 
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Figure 3. Geometry of partial eclipse phases, showing the angles a and /3 used in computing the visible fraction rj of the planetary disc. For illustrative clarity 
only, the planet is shown in front of the star. 



Table 2 

Target pixel positions and fractions of rejected images for each of our 
Spitzer IRAC light curves. The pixel ranges con'espond to the full rang e of 
the Gaussian-filtered positions used for decorrelation in Sect. 14. II 



Target 


Band 


X pixel 


y pixel 


rejected 




^m 


median 


range 


median 


range 


% 


WASP- lb 


3.6 


192.14 


0.15 


152.44 


0.22 


2.1 




4.5 


187.14 


0.15 


94.72 


0.21 


1.8 




5.8 


185.51 




153.57 




3.5 




8.0 


186.26 




94.09 




3.7 


WASP-2b 


3.6 


203.81 


0.06 


149.78 


0.16 


0.3 




4.5 


107.37 


0.05 


155.02 


0.10 


2.6 




5.8 


197.19 




151.04 




3.3 




8.0 


106.22 




153.21 




2.5 



can be seen clearly in the 3.6 fim light cur ve of WASP- lb (top 
left panel of Fig.[Jl. Knut son et alJ (12008!) modelled these sys- 
tematics using a two-dimensional polynomial fit to the instan- 
taneous pointing offset from the centre of the pixel on which 
the stellar image was centred. For observations in which the 
pointing was changed during an observation it was found that 
a quadratic function was needed to remove this intra-pixel ef- 
fect (^nutson et al. 2008; Charbonneau et al. 2008). For ob- 
servations where the pointing position is not changed (like our 
own) some authors find acceptable fits iising j i ist a linear func- 
tion of pixel position jKnutson et al.1 l2009t iMachalek et alj 
I2OO9I: lFressin"eFn]l2(m 

For our analysis we tried both qu adratic and linear dec orre- 
lation functions, similar to those of iKnutson et al.l (l2008h and 
iKnutson et al. (2009). The quadratic function is 

p' = l+cJx-xo)+c_,,,(x-xof+Cy(y-yo)+Cyy(y-yof+c,t (1) 

where (x,y) are the instantaneous coordinates of the stellar 
image, ixo,yo) are the coordinates of the centre of the fidu- 
cial pixel, and t is the mid-time of the exposure. The last 
term allows for a Un ear trend in the light curves as noted by 
IKnutson et al.l(l2009h . The linear decorrelation function is the 
same as Eqn.[Tl except that the constants and c,.,. were set 
to zero. The measured target positions were first smoothed us- 
ing a moving Gaussian filter with a width of 12 observations. 
The median pixel positions andpixel ranges for each of our 
observations are given in Table|2] 

Despite the use of these decorrelation functions, we were 
unable to find acceptable fits when including the steep in- 
crease in flux seen at the beginning of the 3.6 /J,m light curves 



of both planets (Figs. [T| & |2|i. Inspection of the corre- 
sponding backgroun d hght curve s reve aled variations simi- 
lar to that noted by Knutson et al.l (|2009) and attributed to an 
illumination-dependent sensitivity effect si mila r to that seen 
in the 5.8 and 8.0 /im detectors (see Sect. l4.2l i. For the fits 
presented in this paper we decided to exclude the first 30 min 
and 21 min respectively of the 3.6 /im light curves of WASP- 
Ib and WASP-2b. 

4.2. 5.8 and 8.0 jim detrending 

The 5.6 ^m and 8.0 /im (Si:As) detectors suffer from time- 
dependent sensitivity variations that depend on the recent illu- 
mination history of each pixel (Knutson et al. 2007). This can 
be seen clearly in the 8 /im light curves of both WASP- lb and 
W ASP-2b (Figs. [I]&|2ll. Again we followed the methodology 
of lKnutson et al.i t2008t) in fitting a quadratic ramp function in 
logarithmic time, of the form 



p' = l+C,T + C,tT^ 



(2) 



where r = ln(f - fo), f being the mid-time of the exposure and 
to being a fiducial time 30 minutes prior to the first exposure 
of the observing sequence. 

4.3. Secondary eclipse profile 

The secondary eclipse profile is computed assuming the 
planet's day-side hemisphere to have a uniform surface bright- 
ness, and a star/planet flux ratio / = Fp/F^,. 

Given the ratio p = Rp/R^, of the planetary to the stellar 
radius, and a dimensionless separation z in units of the stellar 
radius, the planet is partially eclipsed when l-p < z < l+p. 
At such times, the visible fraction of the planet's disk is given 
by 

/? - cos /3 sin /3 + (a - cos a sin a)/ »^ 
viP^z) = (3) 



where 



and 



cos a = 



l-p^ + z" 



cos/3 : 



2z 



2 2 
-p -z 



2pz 



(4) 



(5) 



using the geometry sketched in Fig. [3] Outside transit, when 
Z > l+p, the visible fraction of the planetary disk is 77 = 1 . 
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Figure 4. Spitzer IRAC light curves covering the secondary eclipses of the exoplanets WASP- lb and WASP-2b. Instrumental effects modelled in our fitting 
process have been removed and the light curves have been scaled to the flux of the star in each band. The light curves have been binned into five hundred phase 
bins per orbital period. 



WASP— lb 3.6/im quadratic decorrelation 




0.45 0.5 0.55 

Orbital phase 

Figures. The 3.6 ^m light curve of WAS P-lb fitted using the quadratic 
decorrelation function described in Sect. 14. II The data have been binned into 
five hundred phase bins per orbital period. Instrumental effects have been 
removed in the bottom panel. 



When z < 1- p, the planet is totally eclipsed and rj = 0. At 
any time, the total observed flux from the star and planet is 
F = F,{l+f7^ip,z)). 

4.4. Parameter fitting 

We solved for the full set of orbital and photometric param- 
eters using the Markov-ch ain Monte-Car lo (MCMC) c ode de- 
scribed by .ColUe r Cam eron et al.l (12007a) and P oUacco et al.l 
(l2008h . modified to incorporate the se condary -eclipse and 
decorrelation models described in Sects. 143144.31 

At each step in the Markov chain, synthetic optical light 
curves and radial-velocity curves are computed using the 
methodology described by PoUacco et al. (2008), using the 
first nine parameters listed in Table[3] This particular set of 
parameters is chose n for thei r mutual near-orthogonahty, a s 
described by Ford (20051) and lCoUier Cameron et al] (l2007al) . 
The secondary eclipse models for th e four IRAC detectors are 
computed as described in Sect. l4.3l multiplying the trial value 



of the planet/star flux ratios / in each bandpass by the planet 
visibility function j] at each time of observation. The normal- 
ized light curve ( 1 + /r/) is then co-multiplied by the appro- 
priate sub-pixel (Eqn.[l]i or ramp (Eqn.|2]i sensitivity model 
with a trial set of model coefficients {cx,Cy,c,} (for linear 
decorrelation in the 3.6 & 4.5 /xm bands), {Cf,c„,Cv,Cyy,c,} 
(for quadratic decorrelation), or {c,,c„} (in the 5.8 & 8.0 /^m 
bands) for each detector in turn. 

The normalized model light curve thus has the form 

p = il+ fi]){ 1 + cAx - xq) + Cyiy - yo) + c,t) (6) 

for linear decorrelation of the 3.6 fjm and 4.5 pm light curves, 

p = {\+fri){\+c_,{x-X())+c^,,{x-XQf+Cy(y-y())+Cyy{y-yQf+c,t) 

(7) 

for quadratic decorrelation of the 3.6 /im and 4.5 /im detec- 
tors, and 

p = {l+flUl+C,T + C„T^) (8) 

for the 5.8 pva and 8.0 pm. light curves. 

The observed fluxes di and the normalized model data 
Pi are orthogonalized by subtracting their respective inverse 
variance-weighted mean values d and p, then computing the 
inverse variance-weighted scale factor 



(di - d)ipi - p)wi/ ^{pi - pfwi. 



(9) 



where the weights w,- = 1 /Var(ii,) are the inverse variances as- 
sociated with the observed fluxes. The logarithmic likelihood 
of obtaining the observed data given the model is quantified 

by _ 

(10) 



X Spitzer " 



Wiidi -d-F^ (pi - p)f . 



This contribution is added to the statistic computed for the 
photo metric and radial- vel ocity data. This is described in de- 
tail bv iPollacco et al.l ( 1200 8). to which the reader is referred. 

At each step in the MCMC calculation, the planet/star flux 
ratios in the four IRAC bands and the nine other model pa- 
rameters {To,P,tTT^F,b,K,ecosuj,esini;j} are each given a 
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Table 3 

Proposal parameter values derived from MCMC parameter fitting to combined optical light curves, radial-velocity curves and Spitzer secondary-eclipse data. 



Parameter 



BJD of primary mid-transit 

Orbital period 

Eclipse duration 

Planet/star area ratio 

Impact parameter 

Stellar mass 

Stellar RV amplitude 



Flux ratio at 3.6 /im 
Flux ratio at 4.5 /im 
Flux ratio at 5.8 /im 
Flux ratio at 8.0 /im 



Symbol 



WASP- lb 

linear decor. quadratic decor. 



To 2453998. 1924 ±0.0002 

P 2.5 19954 ±0.000006 

tr 0.1 550 ±0.0006 

AF 0.0101 ±0.0001 

h 0.064!0U^| 

M, 1.217±d.015 

K 0.11I±0.010 

ecosuj -0.0026 ± 0.0007 

esinw -0.0053 ± 0.0076 

/3 6 0.00184 ±0.00016 

/4 5 0.00217 ±0.00017 

/s 8 0.00274 ± 0.00058 

fg 0.00474 ± 0.00046 



2453998. 1924 ±0.0002 
2.519959 ±0.000006 
0.1550 ±0.0006 
0.0101 ±0.0001 

1.216±0.015 
0.111 ±0.009 
-0.0012 ±0.0007 
-0.0083 ± 0.0073 
0.001 17 ±0.00016 
0.00212 ±0.00021 
0.00282 ± 0.00060 
0.00470 ±0.00046 



WASP-2b Units 
linear decor. 

2454002.2754 ± 0.0002 d 
2.152225 ±0.000003 d 
0.0752 ± 0.0009 d 
0.0178 ±0.0004 

Q 731+0.018 

-0.022 

0.861 ±0.022 M0 
0.154 ±0.004 kms-' 

-0.0013 ±0.0009 
-0.048 ±0.021 

0.00083 ±0.00035 

0.00169 ±0.00017 

0.00192 ±0.00077 

0.00285 ± 0.00059 



random Gaussian perturbation, with a "jump length" of or- 
der the uncertainty in the fitting parameter concerned. The 
Spitzer secondary-eclipse data are then divided by the scaled 
model secondary-eclipse profile, and the parameters of the 
decorrelation function are fitte d by linear least-s quares, using 
singular-value decomposition ('Press et alJll 993h. The proce- 
dure is the same as that used by Gillon et al.l (120101) . The 
model data are computed and fitted to the observations, and 
the global statistic is computed. A decision is then made 
to either accept or reject the proposed parameter set accord- 
ing to the Metropolis-Hastings algorithm: if has decreased, 
the step is accepted unconditionally. If x^ has increased by an 
amount Ax^, the proposal may be accepted with probability 
exp(-AxV2). 

For the first few hundred steps of a typical run, the param- 
eter values evolve in a way that drives x^ down to a global 
minimum value. The algorithm then explores the joint poste- 
rior probability distribution of the parameters in the neigh- 
borhood of th e optim al solution. We use the approach of 
iKnutson et alj (l2008h . declaring the initial burn-in phase is 
complete when x^ exceeds for the first time the median value 
of all previously-accepted values of x^- At this point, we 
rescale the error estimates on the data points in each distinct 
photometric data set, such that the contribution of that data set 
to x^ is approximately equal to the number of observations in 
the set. We then run the chain for a further few hundred steps, 
and compute new jump lengths for the individual fitting pa- 
rameters from the chain variances. After a further short burn- 
in phase, the chain is allowed to continue for a production run 
of 30000 steps. The ensemble of models in this chain defines 
the joint posterior probability distribution for the full set of 
parameters. 

The resulting set of fitting parameters, and the median and 
one-sigma errors of their posterior probability distributions, 
are listed in Table[3] omitting only the rather lengthy list of 
decorrelation coefficients. 

5. RESULTS 

The models best fitting our combined secondary-eclipse, 
primary-transit and radial-velocity data, using linear decorre- 
lation in the 3.6 fim and 4.5 /im bands, are plotted in Figs.[T]& 
|2] The Spitzer data and best-fitting models are also plotted in 
Fig. m where the decorrelation functions have been removed 
for clarity. Fitted parameters are given in Table [3] In common 
with several recent studies, in which the spacecraft pointing 
was not changed during an observation CKnutson et al...2009; 



iMachaleketaTI 120091: iFressin et alJ 120101 ). we find that lin- 
ear decorrelation provides generally a good description of the 
intra-pixel effect in the 3.6 /im and 4.5 /im bands. The ex- 
ception is the 3.6 /im light curve of WASP-lb, where some 
features possibly related to the intra-pixel effect remain. 

We fitted our data also using quadratic de corr elation func- 
tions in the 3.6 /im and 4.5 fim bands (Sect. 14.1b . This made 
no significant difference to the fits of any of our datasets apart 
from the 3.6 /im band of WASP-lb, which is plotted in Fig.|5] 
The fit in this band is improved by the quadratic decorrelation, 
and the eclipse depth is significantly reduced. The best-fit val- 
ues for all parameters from this fit are given in Table |3] 

Our measured secondary-eclipse d epths are compared \y ith 
the planetary atmosphere models of iFortnev et alJ (120081) in 
Fig.|6] For WASP-lb the model predicts a stratospheric tem- 
perature inversion and the two curves represent calculations 
in which the planetary emission is either limited to the day 
side (upper curve) or is uniform over its entire surface (lower 
curve). For the 3.6 /im band we include the eclipse depths 
from both the linear decorrelation (circle) and the quadratic 
decorrelation (square). For WASP-2b, which is much less 
strongly irradiated, the model does not predict a temperature 
inversion and the two curves represent calculations for day- 
side emission (upper) and uniform emission (lower) as before. 
In order to aid comparison with measured eclipse depths we 
have calculated average planet/star flux ratios for each model 
in each IRAC band. These averages have been weighted using 
the IRAC spectral response curves (Hora et al. 2008), which 
are plotted at the bottom of each panel. The model averages 
in each band are represented by horizontal lines. 

The 4.5 /im and 5.8 /im eclipse depths of WASP-lb are con- 
sistent with the model predicting a temperature inversion in 
this system and where emission is limited to the day side of 
the planet. The 5.8 /im band is consistent also with uniform 
emission, but the 4.5 /im band is not. The 3.6 /im eclipse 
depth with linear decorrelation is also consistent with the in- 
version model and day-side emission, but the same data fit- 
ted with quadratic decorrelation is not. Since the quadratic 
decorrelation provides a better fit to the data (Fig.|5]l, this sug- 
gests that the model is over-predicting the planetary flux in 
this band. Alternatively, it may be that the additional degrees 
of freedom in the quadratic decorrelation results in an under- 
estimate of the eclipse depth in this band. 

The 8.0 /im eclipse depth of WASP-lb lies significantly 
above the model prediction (3.3<t). The inversion model pre- 
dicts emission from water in this band, so it may be that the 



8 



Wheatley et al. 




Wavelength (/im) 




Wavelength (/J.m) 



Figure 6. Secondary eclipse depths of the exoplanets WASP- lb and WASP-2b overlaid on the planetary atmosphere models of IFortnev et alj 420081) . which 
are plotted as fraction of the stellar flux. In both cases the upper curves represent models in which emission is only from the day-side of the planet, and the 
lower curves represent uniform emission from the entire exoplanet surface. For the 3.6 ^m band of WASP-lb we have included eclipse depths from both linear 
decorrelation (circle) and quadratic decorrelation (square). The horizontal lines on each model represent averages in the IRAC bands weighted by the spectral 
responses of the camera, which ai'e indicated with solid curves at the bottom of the plot. 



model is under-predicting the strength of the temperature in- 
version in this system. 

In WASP-2b, the ecHpse depths in all four bands are con- 
sistent with the model predicting no temperature inversion 
and with emission only from the day side of the planet. The 
3.6 /im and 5.8 /im eclipses are consistent also with the uni- 
form emission model, but the 8 fim eclipse is only marginally 
consistent with this model and the 4.5 /im eclipse is inconsis- 
tent (3.6(t). 

As further illustration of the importance of the detected 
temperature inversion in WASP-lb, we have also compared 
the measured eclipse depths with a model calculation in which 
the temperature inversion has been artificially suppressed (by 
setting the TiOA'^0 abundances to zero in the model). The re- 
sults are plotted in Fig. |7] and show that the measured eclipse 
depths strongly favor the model with a temperature inversion. 

Table |3] includes updated system parameters for both plan- 
ets from our simultaneous analysis of the secondary-eclipse, 
primary-transit and radial-velocity data. Of particular interest 
is ecosoj, which is constrained by the timing of the secondary 
eclipses. A significant eccentricity is found from the lin- 
ear decorrelation fits to the WASP-lb light curves, however, 
this result is not supported by the quadratic decorrelation fits, 
where the deviation from zero drops to less than two-sigma. 
This fit places a tight three-sigma upper limit of jecoswl < 
0.0033, suggesting that the inflated radius of WASP-lb is un- 
likely to be due to tidal heating. We place a similar upper limit 
on the eccentricity of WASP-2b of |ecosw| < 0.0040. 

6. DISCUSSION 

Overall our measured secondary eclipse depths of WASP- 
lb and WASP-2b in the four IRAC bands a re in remark- 
ably good agreement with the predictions of the' Fortnev et al.l 
(j2008) model (Fig. |6j note that the model has not been fitted 
to the data). This model predicts an atmospheric temperature 
inversion in WASP-lb (class pM) and no inversion in WASP- 
2b (class pL). In WASP-2b the agreement between the model 
and data is excellent in all four bands. In WASP-lb the depth 
of the 8 fim eclipse indicates that the strength of the tempera- 
ture inversion may be under-estimated in the model. In both 
planets the 4.5 /xm eclipse is the best defined and also the one 
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Figure 7. Secondary eclipse depths of the exoplanet WASP-lb overlaid on 
planetary atmosphere models in which the temperature inversion has been 
artificially suppressed by setting the TiOA'O abundance to zero. As in Fig- 
ure|6] the upper curve represents a model with emission only from the day 
side of the planet, and the lower curve represents a model in which heat is 
transported efficiently around the planet. For the 3.6 fim band of WASP-lb 
we have included eclipse depths from both linear decorrelation (circle) and 
quadratic decorrelation (square). The horizontal lines show the weighted av- 
erages of the models in the four IRAC bands, with the IRAC band passes 
indicated by the solid curves at the bottom of the figure. 



least affected by instrumental effects, and in both cases the 
eclipse depth in this band is consistent only with the day-side 
emission model. We take this as strong evidence that redistri- 
bution of incident energy is inefficient in both planets. 

Our detection of a temperature inversion in WASP-lb but 
not in WASP-2b is consistent with expectations based on 
their different irradiation levels (2.5 x 10^ergs"'cm"^ and 
0.9 X 10''ergs"'cm"^ respectively: IFortnev etal] 120081) . It 
is also generally consistent with Spitzer measurements of 
other planets. Of the systems studied in all four IRAC 
bands, HD 189733b, XO-2b and HAT-P-lb have irradiation 
levels at or below the expected transition of around 0.8 x 
10"^ ergs"' cm"', and do not show strong temperature inver- 
sions, although XO-2b and HAT-P-lb do exhibit evide nce for 
weak inversions CCharbonneau et al., ,2008; .Machalek et alj 
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I2OO9I: iTodorov et alJl20Toh . HD209458b, TrES-4 and TrES- 
2 have irradiation levels at or above the expected transition 
and all exhibit evidence f or stronger temperature inversions 
(iKnutsonetalJ 120081 12009: O'Donovanet al. 2010). How- 
ever exceptions to this rule have also been found. XO-lb ex- 
hibits evidence for a temperatu re inversion despite irradiatio n 
at only 0.5 x lO'^ergs"' cm'^ (iMachalek et all l2008l l2009h . 
while TrES-3 has not been found to do so, despite high ir- 
radiation at 1.6 X 10^ ergs"' cm"^ (iFressin et al.ll20ldl) . 

To some extent the observational picture is confused by the 
different models applied by different authors, and the different 
criteria appUed f or the detection of a temperature inversion. 
iGillon et aD (1201 0) plotted a color-color diagram for planets 
studied in all four IRAC bands (their Fig. 10) and found that 
the planets studied to date cover a wide range of color space, 
and do not fall neatly into two distinct groups. In this paper we 
have presented a uniform analysis of secondary eclipse data 
from two planets for the first time. It may be that the obser- 
vational picture will become more clear once further uniform 
analyses of multiple systems are carried out. 

In order to search for a pattern in published secondary 
eclipse depths, we have plotted the ratio of depths in the 
8.0 /im and 4.5 /im bands against irradiation in Fig. |8] We 
chose to compare the ratio of these two bands since they (1) 
are available for the largest number of planets, (2) are sensi- 
tive to the strength of water emission/absorption at 8 fim, (3) 
are less affected by instrumental systematics than the 3.6 /im 
band. Our figure does not show a sharp transition at an ir- 
radiation level of around 0.8 x 10*^ erg s~' cm~^ as m ight be 
expected from the predictions of iFortnev et al.l (120081) . but it 
does show a gradual decrease of the flux ratio through this 
range. However, we do find tentative evidence for a sharp 
transition at a higher irradiation level of 2 x 1 0^ ergs"' cm~^. 
The s ystems beyond this transition are TrES-4 (iKnutson et al.l 
l2009h and WASP- lb (this paper), both showing remarkably 
strong emission at 8 /im, perhaps related to their high irradia- 
tion levels. 

This transition could be due to a change in the chemistry 
and opacity of the atmospheres at that level of incident flux. 
One possibility is that TiO gas is aloft at millibar pressures 
only at higher irradiation levels (and hotter atmospheres) and 
at lower irradiation (< 2 x 10''ergs~' cm~^) another absorber 
is causing the observed inversions, such as a sulfur photo- 
chemical product dZahnle et al. 2009). Both may operate over 
some range. It now appears likely that opacity due to gaseous 
TiO and VO is not the entire story. In particular, the atmo- 
sphere of XO-lb appears to be so cold that Ti should be se- 
questered into solid condensates at any location in the atmo- 
sphere. However, even at slightly warmer temperature in the 
upper atmosphere, there is the still the important issue of the 
temperature of the deep atmosphere. 

Much attention has been paid to how cold traps may af- 
fect the abundance of TiO gas at millibar pressures. If Ti is 
trapped in solid condensates deep in the atmospheres at '--^ 0. 1 
to 1 kbar, then it would require quite strong vertical mix- 
ing to enable gaseous TiO to exist in th e upper a tmosphere. 
This has been discussed by many authors (iHubeny et al..2003; 
IFortnev et aP 120061 2008; Burrows et al. 2008) and calcu- 
lated in som e detail by Showman e t al.i (i2009 ) and especially 
ISpiegel eTa l. (2009). For instance, there is good agreement 
that if TiO causes the temperature invers ion for H D209458b, 
it must be mixed above the cold trap (IFortnev et al. 2008t 
ISpiegel et al.ll2009l) . This is not a priori impossible. In Nep- 
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Figure 8. The ratio of secondary-eclipse depths at 8 /im and 4.5 fim for all 
planets published to date (in planet/star contrast units), plotted as a function 
of irradiation by the parent star. The most highly irradiated planets, including 
WASP- lb, appear to be relatively brighter at 8 fim. 

tune, abundant gaseous methane is found in the stratosphere, 
despite a cold trap that should condense methane in the tro- 
posphere (iFletcher et al.ll20l"0i and references therein). How- 
ever, it is not yet clear whether this mixing through the cold 
trap happens in hot Jupiters._ 

Interestingly, Fortney et al.l (|2006() investigated the deep at- 
mosphere of HD 149026b, which is at an irradiation level of 
2 X 10''ergs"' cm~^ and found that its atmospheric tempera- 
tures sat on the cold trap dividing line (see their Fig. 2), which 
is a tantalizing suggestion that the break seen in Fig.[8]could 
be due to TiO. At higher irradiation levels there is no cold 
trap issue, although the exact location of this boundary is sen- 
sitive to the planet surface gravity and location of the radia- 
tive/convective boundary. We note that Spiegel et al., largely 
because they assumed colder interior temperatures, found that 
this division happened at higher irradiation levels. Additional 
data for the hottest objects will help to support or refute the 
trends shown in Fig. [8] 

7. CONCLUSIONS 

The measured secondary-eclipse depths of WASP- lb and 
WASP-2b presented here indicate a strong temperature inver- 
sion in the atmosphere of WASP- lb, but no temperature inver- 
sion in WASP-2b. This difference is likely to be related to the 
much higher level of irradiation of the atmosphere of WASP- 
Ib. The eclipse depths of both planets also favor models in 
which incident energy is not redistributed efficiently from the 
day side to the night side of the planet. We do not find sig- 
nificant eccentricity in the orbits of either planet, suggesting 
that the inflated radius of WASP- lb is unlikely to have arisen 
through tidal heating. Finally, we find evidence for a sharp 
transition in the properties of planetary atmospheres with irra- 
diation levels above 2 x 10^ ergs"' cm~^, and we suggest that 
this transition might be due to the presence of TiO in the upper 
atmospheres of the most strongly irradiated hot Jupiters. 

Facilities: Spitzer. 
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